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I. INTRODUCTION 



QCD Corrections to Flavor Changing Neutral Coupling Mediated Rare Top Quark 

Decays 
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' J. Stefan Institute, Jamova 39, P. O. Box 3000, 1001 Ljubljana, Slovenia 
'^Department of Physics, University of Ljubljana, Jadranska 19, 1000 Ljubljana, Slovenia 

(Dated: August 5, 2010) 

Recently we have presented an analysis of flavor changing neutral coupling mediated radiative 

top quark decays at next-to-leading order in QCD. In the present paper we provide the details 

of the calculation of QCD corrections to i — >■ 57 and t — >■ qZ decays within the effective theory 

approach including operator mixing. In particular, we calculate virtual matrix element corrections 

and the corresponding bremsstrahlung contributions. In the case of i — > 57 we study the effects 

of kinematic cuts on the extracted branching ratios. Analytical formulae are given at all stages of 

the calculation. We find that the i — > 57 decay can be used to probe also the effective operators 

^^ mediating t ^ qg processes, since these can naturally contribute 10% or more to the radiative decay, 

^~~^ given typical experimental cuts on the decay kinematics at hadron colliders. Conversely, we argue 

^—^ that any positive experimental signal of the i — >■ gg process would indicate a natural lower bound 

^^ on i — >■ 57 decay rate. 

^ PACS numbers: 12.15.Mm,12.38.Bx,14.65.Ha 

Oh 

I The standard model (SM) predicts highly suppressed flavor changing neutral current (FCNC) processes of the top 

M-i quark (i — >■ qV, V = Z,^,g,q = c, u) while new physics beyond the SM (NP) in many cases lifts this suppression (for a 
f-H recent review c.f. [1]). It has been pointed out recently, that top quark FCNC phenomenology is crucial in constraining 

I— —I a wide class of NP scenarios, where new flavor structures are present but can be aligned with the SM Yukawas in the 
down sector [2H5]. Top quark FCNCs can be probed both in production and in decays. Presently the most stringent 

%^ bound on the Br(i — > qZ) comes from a search performed by the CDF collaboration Br(i — > qZ) < 3.7% at 95% 
C.L. [B^. The photonic decay is presently mostly constrained by the ZEUS collaboration Br(i — ^ qj) < 0.59% at 95% 
C.L. |7j. On the other hand the most stringent present limit on Br(i -^ qg) comes from the single top production 

ly-^ total cross-section measurement of CDF and yields Br(i -^ ug) < 0.039% and Br(i — >• eg) < 0.57% at 95% C.L. [S]. 

^Nl The LHC will be producing about 80,000 ii events per day at the luminosity L = lO^'^cm^^s"'^ and will be able to 

P I access rare top decay branching ratios at the 10"*^ level with lOfb^^ [9^. 

(^ Recently [inilll] the t -> qV decays mediated by effective FCNC couplings have been investigated at next-to-leading 

f^ order (NLO) in QCD. In ^Jj it was found that t -^ qg receives almost 20% enhancement while corrections to the 

T-H t — > C7, Z branching ratios are much smaller. Contributions of additional operators in i — > qZ and the effects of 

operator mixing induced by QCD corrections have been identified in [TD] . In the case of t —> q-f decay in particular 
the QCD corrections generate a nontrivial photon spectrum and the correct process under study is actually t — >■ qjg. 
Experimental signal selection for this mode is usually based on kinematical cuts, significantly affecting the extracted 
bounds on the effective FCNC couplings. The main implications of this observation were presented in [10 , while the 
details of the underlying calculation and full analytic results are the subject of the present paper. 

This paper is structured as follows. We begin by giving the set of FCNC operators considered and explain our 
notation. Next we conduct a full study of QCD corrections to t — ^ qZ, 7 decays. Contributions from gluonic dipole 
operators are also taken into account. We present results for the virtual matrix element corrections as well as the 
corresponding bremsstrahlung rates. In the i — )■ (77 channel we also study the relevance of kinematical cuts on the 
photon energy and the angle between the photon and the jet stemming from the final state quark. We present our 
results in analytical form and also give numerical values to estimate the significance of NLO contributions. 
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II. FRAMEWORK 



In writing the effective top FCNC Lagrangian we rely on the notation of ref. [U [T2]. Hcrmitian conjugate and 
chirality flipped operators are implicitly contained in the Lagrangian and contributing to the relevant decay modes 
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To explain the notation, operators considered are 



Of.R = 9zZ^ 



^LR,RL — '-^tJ.'J 



(lL,Rl^tLM^ 
qL,R(j'"'tR^L 



^LR,RL — gzZ^j,^ 



qL,Rcr^''tR^L 
QL.R'^'^^TatR^L 
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where qr l = (1 ± 75)9/2, cr^iy = i['y^,ji,]/2 and gz = 2e/ s'm29\Y- Here 9w stands for the Weinberg angle, while 
e = V4^ and g, = V^i^- Furthermore V{A, Z)^, = d^V, - d,V^ and CJ^^ = d^G^ - d^C^ + gfabcGlGl where 
and T" and fabc are the SU{'i) color group generators and structure constants respectively. Finally v = 246 GeV is 
the electroweak condensate and A is the effective scale of NP. In the remainder of the paper, since there is no mixing 
between chirality flipped operators we shorten the notation, setting a and h to stand for ai,, hLR or ur, Brl. 

Note that in principle, additional, four-fermion operators might be induced at the high scale which will also give 
contributions to t — )■ qV processes, however these are necessarily a^ suppressed. On the other hand, such contributions 
can be more directly constrained via e.g. single top production measurements and we neglect their effects in the present 
study. We also assume the effective a, b couplings are defined near the top quark mass scale at which we evaluate 
virtual matrix element corrections and a^. A translation to a higher scale matching is governed by the anomalous 
dimensions of the effective operators and can be performed consistently using RGE methods. The procedure and its 
implications has been discussed in ref. [10) and we do not repeat it here. In our calculation we neglect the mass of 
the final state (c, u) quark and regulate UV as well as IR divergences by working in d = 4 + e dimensions. 

III. t -^ qZ DECAY 

The total t -^ qZ decay width, mediated by operators in Lagrangian (IT]) and including leading QCD corrections 
can be written in the following form 
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l^^prf + |6SpFf + 2Re{6^*&9}Ffg - 2Im{6^*6S}ffg 



(3) 



It includes tree level contributions, one loop virtual corrections and gluon bremsstrahlung processes (actually t — >■ qZg) 
which we will analyze separately in the following two subsections. In the end of the section we present numerical 
results of the combined effects. 



A. Tree level expressions 

At the tree level we only have F^, F^ and F^j^ contributions, which we write in 4 + e dimensions as 

Ff(°) = lmr^g|(l-r^mi+2)(l--^r2^(l + (2 + ^W)> 
rf'"^ = lim^.g|(l-rzmi + ^)(l-rzr2(2 + e + rz), 
rfr = hm^.9|(l-r^mi + |)(l-r;,r(3 + e). 



(4) 



where rz = to^/toj . 



B. Virtual corrections 



At the one loop level, t — > qZ, 7 decay rates receive contributions from Feynman diagrams in Figure [Tl Besides 
the one loop gluon corrections to the Z, 7 FCNC operator matrix elements, we also include contributions from the 
gluonic FCNC dipole operators. The later first appear at 0{as) with the emission of the Z,j proceeding through SM 
couplings. On the other hand, one particle reducible diagrams with gluon corrections attached to the external legs 




Figure 1: One loop virtual corrections to i — >■ cZ, 7 decay. Crosses mark additional points from which the Z, 7 boson can be 
emitted. 

are not presented in Figure 1 They are taken into account by proper quark field renormalization q — > ^TZ^q. Since 
the final state light quark is considered to be massless, the corresponding field renormalization differs from the one of 
the initial top quark. Using the on-shell renormalization conditions we get 
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where fi is the renormalization scale parameter and C'f = 4/3. The resulting t — >■ qZ decay amplitude including one 
loop virtual effects can be written in terms of six form factors 
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The form factors are 
F^ = C, 
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We have defined auxifiary functions /i and /2 for shorter notation 
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F^ and i^^^ include the quark field renormalization. To rid all the UV divergences, operator renormalization is 
necessary. This leads to the appearance of counter terms 6^ and S^ which are renormalized in the appropriate 
matching condition in the UV. The two operator renormalization counter terms are evaluated in MS scheme 



r,Z 



yZ 



47r 



Cpsi 



CpSbg ■ 



5f 



< 



(- 

^euv 



+ 7 - log(47r)) 



+ 7 - log(47r) 



(10) 
(11) 



On the other hand, C^ = {rrit/ ^YT{1 — e/2)/(47r)'" is an IR renormalization factor multiplying divergent form 
factors. Finally, SM Z boson couplings to fermions are defined as w = T3 — 2sin6'vvQ and a — T3, where for the 
up-type quarks T^ — 1/2 and Q — 2/3. They appear only in form factors generated by the gluonic dipole operator. 

We were able to crosscheck our expressions for form factors with those found in the literature. Namely, (7|) agree 
with the corresponding expressions given in ref. [13] for the B — >■ Xsl^l~ decay mediated by a virtual photon 
after taking into account that the dipole operator in [T3] includes a mass parameter which necessitates additional 
mass renormalization. To some extent we were also able to crosscheck the gluon operator induced form factors (|8|. 
Namely, we find numerical agreement of the form factor vector component (i.e. setting d = 0) with the corresponding 
expressions given in ref. [14]. The crosscheck is of course only possible in the vector part, since the SM photon 
coupling appearing in ^14] has no axial component. 

After UV renormalization we are still left with IR divergences. These will be canceled at the level of the decay width 
by the IR divergences associated with the bremsstrahlung process t — > qZg. In order to demonstrate this cancellation 
explicitly, we write down the contribution of the tree level and virtual corrections to decay widths defined in eq. (|3| . 
The IR divergent parts are 
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Remaining contributions are induced by the gluonic operator and are IR finite 
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C. Bremsstrahlung Contributions 

The relevant Feynman diagrams contributing to f — > qgZ, 7 bremsstrahlung processes are given in Figure [2J At the 
level of the decay width these diagrams give contributions of the same order in a^ as the one-loop virtual corrections 
presented in the last subsection. Soft and collinear IR divergences emerge in the phase space integration. Below 




Figure 2: Feynman diagrams for t — >■ cgZ, 7 bremsstrahlung process. Crosses mark additional points form which the gluon or 
Z, 7 boson can be emitted. 

we give analytical formulae for the bremsstrahlung contributions to the decay widths in eq. (|3|. There arc three IR 
divergent contributions 
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Summing these decay widths with the one-loop virtually corrected t —> qZ rates into F 
X = a, 6, ab, we obtain combined IR finite expressions 
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Again we were able to crosscheck our results (151 with the corresponding calculation done for a virtual photon 
contributing to the B — >■ Xsl'^l~ spectrum [121 . After taking into account the different dipole operator renormalization 
condition in [15j (including mass renormalization) we find complete agreement with their results. T^ was also cross- 
checked with the corresponding calculation of the t — >■ Wh decay width at NLO in QCD 16^. Finally, we have 
compared our F^ expression with the results given by Zhang et al. in ref. |11| . In the limit rz — >■ our results agree 
with those given in |TT] , but we find disagreement in the rz dependence. After our first publication of these results 
in |10j . we were made aware of a new paper in preparation by the same authors, which has now been published [17] 
and therein a corrected result for T\ is given that coincides with ours. 

The remaining bremsstrahlung contributions are induced by the gluon dipole operator and are IR finite 
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Numerical agreement was found when comparing the vector parts of our resuhs given in eq. ( 16 1 with the corresponding 
results of ref. 1151 . 



D. Numerical analysis 

In this section we present some numerical values to estimate the significance of QCD corrections. In particular we 
parametrize the decay width given in eq. (|3| as 
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Here a;^ stand for the tree-level contributions, while j/i , yi denote the corresponding QCD corrections. Numerical 
values of the coefficients are given in Table [ij We see that corrections due to the gluon dipole operator are an order 



of magnitude smaller (except i/g , which is even more suppressed) than corrections to the Z operators themselves and 
have opposite sign. 
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Table I: Numerical values of coefScient functions at the renormalization scale equal to the top quark mass corresponding to the 
inputs m,t = 172.3 GeV , mz = 91.2 GeV , sin^ Qw = 0.231 . 

Next we investigate the relative change of the decay rates and branching ratios when going from the leading order 
to next to leading order in QCD, where branching ratio for the top quark is defined to be normalized to its main decay 
channel Br(t — >■ qZ,j) = r{t — )■ qZ,^)/r{t — >■ bW). The results are presented in Table [TT| We see that the change in 
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Table II: Numerical values of p^^^'-'/r^" and Br"'"'-'(t -5> qZ)/'Br^^{t ->■ qZ) for certain values and relations between Wilson 
coefficients. Results are obtained using analytical formulae given above setting fj. = rrit. Main decay channel decay width at 
NLO in QCD is given in ref pjj. Additional SM inputs used are niw = 80.4 GeV and q^ = 0.107. 



the decay width is of the order 10%. There is a severe cancellation between the QCD corrections to r(t — >■ cZ) and the 
main decay channel r(i — )■ hW). This cancellation causes the change of the branching ratio to be only at the per-mile 
level when h^ is set to zero. In the case when only operators Of ^ are considered this cancellation is anticipated since 
the NLO correction to Ff is of the same form as the correction to the rate of the main decay channel. If we treat 
b quarks as massless, exact cancellation is avoided only due to the difference in the masses of Z and W bosons. It 
is more surprising that similar cancellation is obtained also when only the dipole Z operator is considered. However, 



setting &s 
percent. 



or b^ = b , the impact of QCD corrections is increased by an order of magnitude and reaches a few 



IV. t -^ q-j DECAY 

Compared to eq. ([3]), the expression for the total t ^f q^/ decay width gets simplified, since due to gauge invariance 
we only need to consider dipole operators 
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We follow a similar pattern of analysis as in the Z case, with the added complication that the bremsstrahlung 
contributions have to be treated in greater detail since the non-trivial photon spectrum has important implications 
for the experimental detection of this decay channel. We start with the tree level contribution, which we write in 4-|-e 
dimensions as 



r 



^(°) = limm,a(l + ^)r(l + ^). 



(19) 



A. Virtual corrections 

Feynman diagrams yielding virtual corrections to the t —>■ q^ decay rate are presented in Figure [l] We are dealing 
with one less diagram, since we only have the photonic dipole operator. Quark field renormalization remains the same 
and the decay amplitude can be written in terms of two form factors 



Af. 



>97 



V 
A2 



b^l+'^Cpr 



^'t.'^^K. 



\^LR,RL/ ' 



(20) 



g 0.6 






1 0.4 




0.0 0.2 0.4 0.6 0.8 1.0 

l-2Ey/mt 



Figure 3: The t — >■ qyg Dalitz plot. Contours of constant photon and gluon infrared and collinear divergent contributions 
are drawn in red (dot-dashed) and blue (dashed) lines respectively. The collinear divergencies appear at the horizontal and 
vertical boundaries of the phase-space, while the IR divergencies sit in the top and right corners. The cuts on the photon 
energy correspond to vertical lines, the cuts on the gluonic jet energy to horizontal lines. Full green lines correspond to cuts 
on the jet veto cone around the photon. 



which read 
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(22) 



Operator renormalization counter terms induced by the UV matching procedure in the MS scheme read 
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Finally, tree level and one-loop virtual correction contributions to the t — > 97 decay width are 
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(23) 



(24a) 

(24b) 
(24c) 



As in the t — > qZ case we are left with IR divergences which have to be canceled by the corresponding bremsstrahlung 
contributions. 



B. Bremsstrahlung corrections 



The t — )■ 5717 decay process involves three (one almost) massless particles in the final state. Virtual matrix 
element corrections contribute only at the soft gluon endpoint (Eg = 0) and result in non- vanishing b'^b^ interference 
contributions. They involve IR divergencies which are in term canceled by the real gluon emission contributions. 
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These also produce non- vanishing |6^p contributions, and create a non-trivial photon spectrum involving both soft 
and coUinear divergencies. The later appear whenever a photon or a gluon is emitted coUinear to the light quark 
jet. An analogous situation is encountered in the B -^ X^'y decay measured at the i3-factories. However, there 
the photon energy in the B meson frame can be reconstructed and a hard cut {E^^^) on it removes the soft photon 
divergence. The cut also ensures that the B -^ Xgg process contributing at the end-point E^ — is suppressed. On 
the other hand, in present calculations the coUinear divergencies are simply regulated by a non-zero strange quark 
mass, resulting in a moderate log(ms/?nf,) contributions to the rate. The situation at the Tevatron and the LHC 
is considerably different. The initial top quark boost is not known and the reconstruction of the decay is based on 
triggering on isolated hard photons with a very loose cut on the photon energy (a typical value being E^ > 10 GeV in 
the lab frame [TSj). Isolation criteria are usually specified in terms of a jet veto cone AR = y/Arf + A^^ where Arj 
is the difference in pseudorapidity and A0 the difference in azimuthal angle between the photon and nearest charged 
track. Typical values are AR > (0.2 — 0.4) [15]. We model the non-trivial cut in the top quark frame with a cut on 
the projection of the photon direction onto the direction of any of the two jets {Sr = 1 — p^ • pj/E^Ej). The effects 
of the different cuts on the decay Dalitz plot are shown in Figure [3) Since at this order there are no photon coUinear 
divergencies associated with the gluon jet, the x cut around the gluon jet has a numerically negligible effect on the 
rate. On the other hand the corresponding cut on the charm jet - photon separation does not completely remove the 
divergencies in the spectrum. However, they become integrable. The combined effect is that the contribution due to 
the gluonic dipole operator can be enhanced compared to the case oi B —>■ Xgj. Below we give the full analytical 
formulae for the t — ^ 575 decay rate including the effects of kinematical cuts on variables x = 6r and y = 2i?™*/r7it. 
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It is easy to verify the cancellation of IR divergences upon summation of V 
(25a I respectively. 
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Figure 4: Relative size of a^ corrections to the Br(i — >■ ^7) at representative ranges of Sr and Ej 
correction values are plotted for 6® = (gray, dotted), b® — b'' (red) and 6* — —6^ (blue, dashed). 
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Figure 5: The ratio of radiative rates V(t — > q'y)/T{t — >■ qg) versus the absolute ratio of the relevant effective FCNC couplings 
\b~' /h^\. Two representative choices of experimental kinematic cuts are shown. The shaded bands represent the possible spread 
due to the unknown relative phase between 6^ and b^ couplings, while the lines correspond to maximal positive (full, dotted) 
and negative (dashed, dot-dashed) interference b^b^ . See text for details. 



C. Numerical analysis 



In Figureplwe show the b^ induced correction to the tree-level Br(< — > 97) for representative ranges of 5r and E"*. 
We observe,^hat the contribution of h^ can be of the order of 10 — 15% of the total measured rate, depending on 
the relative sizes and phases of 0£'^ ^^ ^'^'^ '^'^ ^^^ particular experimental cuts employed. Consequently, a bound 
on Br(i — > (77) can, depending on the experimental cuts, probe both 6^''*' couplings. In order to illustrate our point, 
we plot the ratio of radiative rates r(i — > q^)/T{t — ?► qg), both computed at NLO in QCD versus the ratio of the 
relevant effective FCNC dipole couplings \b'' /b'^\ in figure l5J We show the correlation for two representative choices 
of experimental kinematic cuts for the i — > (77 decay. The vertical spread of the bands is due to the variation of the 
relative phase between b'' and b^ couplings. We also display the two interesting limits where the b'^b^ interference 
is maximal positive (zero relative phase) and negative (relative phase tt). We see that apart from the narrow region 
around \W /b^\ '^ 0.2, where the two contributions may be fine-tuned and conspire to diminish the total t — > C7 
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rate, the two radiative rates are well correlated. In particular, depending on the kinematical cuts employed, there 
is a natural lower bound on ratio of decay rates, valid outside of the fine-tuned region. Finally, for 1 6^/6^1 > 0.6 
the correlation becomes practically insensitive to the particular experimental cuts employed and also the unknown 
relative phase between IP and h^ couplings. 

V. CONCLUSIONS 

In summary, QCD corrections to FCNC coupling mediated rare top decays can induce sizable mixing of the relevant 
operators, both through their renormalization scale running 110] as well in the form of finite matrix element corrections. 
These effects are found to be relatively small for the t ^^ cZ decay. On the other hand the accurate interpretation of 
experimental bounds on radiative top processes in terms of effective FCNC operators requires the knowledge of the 
experimental cuts involved and can be used to probe Cf,^. ^^ contributions indirectly. 

Finally we note that additional information on the underlying NP contributions is provided by kinematical distri- 
butions of i — >■ cZ, 7 common final states such as charged di-leptons [20] as well as in the case oit ^ eg through single 
top production cross-section ^21ti^ . Combined, these observables could facilitate the reconstruction of prospective 
NP models in case a positive experimental signal of FCNC top quark processes would emerge in the future. 

Acknowledgments 

J. F. K. would like to thank Mikolaj Misiak and Gino Isidori for useful discussions and the Galilo Galilei Institute 
for Theoretical Physics for the hospitality and the INFN for partial support during the completion of this work. 
This work is supported in part by the European Commission RTN network. Contract No. MRTN-CT-2006-035482 
(FLAVIAnet) and by the Slovenian Research Agency. 



[1 
[2 

[3: 
[4; 
[s: 
[6: 
[7: 

[8 
[9 

[lo; 
[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23 



J. A. Aguilar-Saavedra, Acta Phys. Polon. B35, 2695 (2004), hep-ph/0409342. 

P. J. Fox, Z. Ligeti, M. Papucci, G. Perez, and M. D. Schwartz, Phys. Rev. D78, 054008 (2008), 0704.1482. 

O. Gedalia, L. Mannelli, and G. Perez (2010), 1002.0778. 

O. Gedalia, L. Mannelli, and G. Perez (2010), 1003.3869. 

A. Datta and M. Duraisamy, Phys. Rev. D81, 074008 (2010), 0912.4785. 

T. Aaltonen et al. (CDF), Phys. Rev. Lett. 101, 192002 (2008), 0805.2109. 

S. Chekanov et al. (ZEUS), Phys. Lett. B559, 153 (2003), hep-ex/0302010. 

T. Aaltonen et al. (CDF), Phys. Rev. Lett. 102, 151801 (2009), 0812.3400. 

J. Carvalho et al. (ATLAS), Eur. Phys. J. C52, 999 (2007), 0712.1127. 

J. Drobnak, S. Fajfer, and J. F. Kamenik, Phys. Rev. Lett. 104, 252001 (2010), 1004.0620. 

J. J. Zhang et al., Phys. Rev. Lett. 102, 072001 (2009), 0810.3889. 

J. A. Aguilar-Saavedra, Nucl. Phys. B812, 181 (2009), 0811.3842. 

A. Ghinculov, T. Hurth, G. Isidori, and Y. P. Yao, Nucl. Phys. B648, 254 (2003), hep-ph/0208088. 

A. Ghinculov, T. Hurth, G. Isidori, and Y. P. Yao, Nucl. Phys. B685, 351 (2004), hep-ph/0312128. 

H. H. Asatryan, H. M. Asatrian, C. Greub, and M. Walker, Phys. Rev. D66, 034009 (2002), hep-ph/0204341. 

C. S. Li, R. J. Cakes, and T. C. Yuan, Phys. Rev. D43, 3759 (1991). 

J. J. Zhang et al. (2010), 1004.0898. 

G. Aad et al. (The ATLAS) (2009), 0901.0512. 

CMS Note 2007^06 . 

J. Drobnak, S. Fajfer, and J. F. Kamenik, JHEP 03, 077 (2009), 0812.0294. 

E. Malkawi and T. M. P. Tait, Phys. Rev. D54, 5758 (1996), hep-ph/9511337. 

M. Hosch, K. Whisnant, and B. L. Young, Phys. Rev. D56, 5725 (1997), hep-ph/9703450. 

T. Han, M. Hosch, K. Whisnant, B.-L. Young, and X. Zhang, Phys. Rev. D58, 073008 (1998), hep-ph/9806486. 



